Coulomb gap in one-dimensional disordered electronic systems 
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We study a one-dimensional system of spinless electrons in the presence of a long-range Coulomb 
interaction (LRCI) and a random chemical potential at each site. We first present a Tomonaga- 
Luttinger liquid (TLL) description of the system. We use the bosonization technique followed by the 
replica trick to average over the quenched randomness. An expression for the localization length of 
the system is then obtained using the renormalization group method and also a physical argument. 
We then find the density of states for different values of the energy; we get different expressions 
depending on whether the energy is larger than or smaller than the inverse of the localization length. 
We work in the limit of weak disorder where the localization length is very large; at that length 
scale, the LRCI has the effect of reducing the interaction parameter K of the TLL to a value much 
smaller than the noninteracting value of unity. 

PACS number(s): 71.10.Pm, 71.27. +a, 71.20.-b 



I. INTRODUCTION 

In a strongly localized electron system with an un- 
screened long-range Coulomb interaction (LRCI) be- 
tween the electrons, the density of states (DOS) is ex- 
pected to exhibit a gap at the Fermi energy. This is called 
the "Coulomb gap" [1] as it arises due to the LRCI. A 
phenomenological argument [1] suggests that the DOS 
should obey a power-law of the form 

D(w) ~ M d - X , (1) 

where d is the spatial dimensionality of the system, and 
uj is measured from the Fermi energy. We shall consider 
here the one-dimensional case where the DOS is expected 
to show a logarithmic behavior in the extremely localized 
limit [2-5]. 

Following the recent developments in carbon nanotube 
technology [6], there has been an upsurge in the studies 
of one and quasi-one dimensional electron systems. In 
low dimensional systems, the electron-electron interac- 
tions play a dominant role leading to a behavior signifi- 
cantly different from that of conventional Fermi liquids. 
Short-range repulsive interactions between the electrons 
lead to the Tomonaga-Luttinger liquid (TLL) behavior 
[7-9], whereas the LRCI is believed to lead to a Wigner 
crystal [10]. 

It is well known that for a one-dimensional noninteract- 
ing electron system with random disorder, all the states 
are localized due to repeated back-scatterings of the elec- 
trons [11]. However, the DOS of noninteracting electrons 
is finite at the Fermi energy [12]. The presence of LRCI 
together with random impurities modifies the DOS in a 
drastic way and supposedly leads to the Coulomb gap 
behavior as discussed above. In an earlier study, Vojta 
and John [5] studied a one-dimensional electron system, 
with LRCI and randomness, in the extremely localized 



limit where the overlap of electronic wave functions can 
be neglected. They find the form of the DOS to be [5] 

D M ~(l»g)"\ (2) 

where E c is a cut-off energy [3,4]. In a recent pa- 
per, Lee [13] derived the expression of the DOS of a 
one-dimensional system of spinless electrons with LRCI 
and impurities at random positions in the opposite limit 
where the quantum effects (electron hopping) play a dom- 
inant role, and the pinning is weak, i.e., the localization 
length is much larger than the inter-impurity distance. 
Following studies of a pinned Wigner crystal [14,15] and 
using a semiclassical approximation, Lee [13] finds the 
DOS at low energy to be 

D(u) ~ , (3) 

where the exponent 7/ (defined in Eq. (19) below) is de- 
termined by the localization length of the system. This 
power-law behavior is consistent with the existing numer- 
ical studies [16]. 

In this paper, we study a model of spinless electrons 
with LRCI as in Ref. [13] (with quantum effects, i.e., 
with the electron hopping term) in the presence of a ran- 
dom chemical potential at each site. In our approach, 
we find that a TLL picture [9], rather than a semiclas- 
sical approach [13], provides a convenient description of 
the system. The localization length Lq is derived us- 
ing a renormalization group (RG) study of the effective 
disorder- averaged bosonized action [17], and is found to 
be very large in the limit of weak disorder. (We also 
present a simple physical understanding of the expression 
for Lq). Due to the LRCI, the interaction parameter K 
of the TLL is found to be effectively a function of the 
length scale. At the length scale L , the value of K is 
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given by a value much smaller than unity (which is the 
value of K for the noninteracting system). We should 
remark here that the present one-dimensional disordered 
electron system with LRCI appears to be an unique ex- 
ample of a TLL with K « 1. 

Once we make use of a TLL description, the DOS can 
be found in a standard way [9,18]. We will show that for 



> 



vp/Lo, the DOS is given by 



D(u) 



LO 







+ 



- 1 



(4) 



where, unlike a standard TLL, rj itself depends on u as 
discussed below. Eq. (4) resembles the result of Lee [13] 
in the limit r] » 1; however, -q is not taken to be a func- 
tion of lo in Ref. [13]. On the other hand, for u> < vp/L , 
we will argue that the DOS is given by the expression in 
Eq. (2). 

In Sec. II, we will consider a system of spinless elec- 
trons with LRCI and disorder. After bosonizing the 
Hamiltonian, we well use the RG equations to obtain the 
localization length L in terms of the disorder strength. 
We will also present a physical understanding of the ex- 
pression for L . In Sec. Ill, we will obtain expressions 
for the DOS for different ranges of values of u. In Sec. 
IV, we will briefly discuss the case of a screened Coulomb 
interaction. 



II. BOSONIZATION AND THE LOCALIZATION 
LENGTH 

The Hamiltonian of a disordered system of spinless 
electrons with LRCI consists of three parts, 



H — H + He + ^random 



(5) 



The noninteracting part H and the Coulomb interaction 
part He are written in terms of the continuum chiral 
electron fields as follows, 



H a =v F / dx (-iip R d x ipR + itpld x ip L ) , (6) 

J — oo 

i r°° 

H c = - dxdy U(x- y)p(x)p(y) , (7) 

Z J-oo 

where the form of U(x — y) will be specified later. The 
fields iPr(iI>l) are the right-moving (left-moving) electron 
operators; they are related to the lattice electron opera- 
tors by 



c(x) = y/E [e lkFX ^ R (x) + e- ikFX ^ L {x)] 



(8) 



where a is the lattice spacing. The density p = Pr + Pl, 
where pn (pp) is the normal-ordered density of the right 
(left) moving electrons. 



In our model, the random part of the Hamiltonian con- 
sists of two parts in the continuum limit, (i) the forward 
scattering part, Hf, where the scattered electrons remain 
in the vicinity of the same Fermi point, and (ii) the back- 
ward scattering part, Ht, where an electron is scattered 
from — kp to kp, or vice versa. 

Using the rules of bosonization [9,18], we rewrite the 
full Hamiltonian in the bosonic language. The low-energy 
and long-wavelength excitations of the noninteracting 
part can be written as 



/oo 
dx [(d x 9{x)f + {d x ^(x)f] , 

while the Coulomb part is given by 

1 r 00 

i7 c = — / dxdy U(x- y)d x 4>{x)d y 4>{y) . 

^ 7r J — oo 

The forward scattering part is given by 

H f = I dx h(x) ip* R {x)i) R {x) + ipl(x)ip L (x) 
dx h(x)d x (j)(x) . 



(9) 



(10) 



i r 

7T J-, 



(ii) 



The forward scattering is due to the real random field 
h(x), and its effect can be taken care of by a rescaling of 
the bosonic field <j> [17]; henceforth, we shall ignore the 
forward scattering part. The back-scattering part can be 
written as 



H b 



dx 



^{x)^ r {x)^ l {x)+C{x)^ l {x)^r{x) 



i r 

7ra J_, 



dx 



£(a;) e *(2*(*)+2*F*) +hx ; 



(12) 



where h.c. denotes the hermitian conjugate. The back- 
scattering is due to a complex random field £ (x) with the 
probability distribution 



P[$(x)] = exp 



-Df j dxC(xMx) 



(13) 



so that < >= and < £*(x)£(x') >= D s 5(x - x'). 

The Hamiltonians given in Eqs. (9), (10) and (12), 
along with the probability distribution (13) of the ran- 
dom field £, constitute the complete low-energy and long- 
wavelength description of the model in the bosonic lan- 
guage. We now rewrite the Coulomb interaction in Eq. 
(10) in momentum space as 



^ J dxdyU{x - y)d x (f>{x)d y <j)(y) 

^3 J dkk 2 U(k)<j>(k)<j>(-k) , 



(14) 



where U(k) is the Fourier transform of the Coulomb in- 
teraction. We can expand U(k) in powers of k as 
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U(k) = U + U 1 k + U 2 k 2 + 



(15) 



From the unperturbed Hamiltonian in Eq. (9), we find 
the naive scaling dimensions of various quantities to be 



[k] = (Length)" 1 , [0] = (Length) 1 . (16) 



Considering the scaling dimensions of the various terms 
in Eq. (15), we conclude that only the constant term is 
marginal under renormalization, whereas the terms in- 
volving powers of k are irrelevant. Hence, the important 
contribution arising from the Coulomb interaction term 
in Eq. (14) can be written as 



H c = ^ J dkk 2 j>(k)j>(-k) 
= S 2 J dx ( d ^( x )f ■ 



The constant Uq is given by 



U = 



L/2 



dx 



L/2 Vx 2 + d 2 



2q 2 In - 

a 



(17) 



(18) 



where q is the charge of the electron, and d is the width 
of the wire which is of the order of the lattice spacing a; 
it ensures convergence of the integral in Eq. (18) near 
x = 0. The meaning of the cut-off length L will be dis- 
cussed below. (We have assumed that L » a,d; this 
will be justified later for various cases of interest). We 
now define a dimcnsionlcss quantity 



2q 2 L 
i] = In — 

ttvf a 



(19) 



Then the Hamiltonian can be written in the form H = 
H + H c + H b , where 



Hq + He 



VF 

2ir 



dx[(d x e{x)f + {i + 7 1 ){d x <t>{x)Y 



(20) 



Eq. (20) can be brought into a standard form by defining 
1 



K = 



We then arrive at the expression 



and u — vpyjl + rj . 



(21) 



H a + He 



u 
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dx 



K{d x e{x)f 



K 



{d x cj>{x)Y 



(22) 



In the imaginary time representation, the quadratic 
and random parts of the action are respectively given by 



So 



Sc 
u 

2^ 
1 

ira 



dxdr 
dxdr 



K 



(d x <P(x,r)) 2 



^ x y{24>{x,r)+2k F x) + hc 



(23) 



The randomness in the back-scattering term is dealt with 
using the standard replica method and averaging over the 
randomness using the distribution in Eq. (13). The fi- 
nal form of the disorder-averaged n-replicated action is 
found to be [17] 



= ^— / dxdr 



K(d T r(x,r)) 2 + -(d x( (> a (x,T)) 2 



D f 



{na) 2 



J dxdrdr' cos(20 a (x,T) -2/(x,t')), 

a.b 



(24) 



where a, b are the replica indices running from 1 to n. 
Eq. (24) is the same as the replicated action of a one- 
dimensional disordered electronic system given in Ref. 
[17]. However, the back-scattering term arising due to 
the spin of the electrons is absent in the present case 
of spinless electrons. The RG equations of the different 
parameters are given by [17] 



dK 

~~dT 

du 



-K 2 D 
2 



uK 



D 



(25) 



where D = (2D ^a) / 'iru 2 , and I — ln(L/a) is the logarithm 
of the length scale. The RG equations given above show 
a quantum phase transition at K = 3/2. For K > 3/2, 
the disorder is irrelevant. For K < 3/2, the disorder 
grows under RG and leads to localization. Thus the sys- 
tem undergoes a zero-temperature localized-delocalized 
transition at K = 3/2. 

The localization length, denoted by Lq, can be ob- 
tained by integrating the RG equations in (25). At the 
microsocopic length scale a, we begin with some values 
of D (much smaller than 1), and u, K given in Eq. 
(21). (The latter requires a knowledge of the length 
scale L through Eq. (19); this will be determined self- 
consistently by the solution of the RG equations). Since 
we are assuming that K « 1, we find (as can be verified 
numerically using Eqs. (25)) that the quantities K and 
u flow very little while D flows from a small number of 
order to a number of order 1. Thus the length scale at 
which the disorder strength becomes of order 1 is given 
by 



Lo 
a 



D f 



l/(3-2if) 



(26) 
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We identify this as the localization length whose physical 
meaning will be discussed below. 

In the limit of weak disorder (D^ — > 0), the form of the 
localization length in Eq. (26) implies, in a self-consistent 
manner, that (i) Lq >> a, (ii) the parameter rj >> 1 
(replacing L by Lq in Eq. (19)), and (iii) the interaction 
parameter K « 1 due to Eq. (21). Hence the local- 
ization length in the K ^ limit assumes the classical 
value 



Lo 
a 



1 



(27) 



[Note that the bosonization approach cannot be used 
in the extremely localized limit (Lq <~ a) studied in Rcfs. 
[3-5]. In that limit, one cannot use a continuum descrip- 
tion, and thus a bosonized description of the system is 
not possible]. 

The significance of the localization length is as follows. 
Although two electrons which are separated by more than 
the distance Lq do interact with each other through the 
Coulomb potential, the overlap of their wave functions is 
exponentially small, and hence their positions are uncor- 
rected. Such interactions will therefore only contribute 
to a uniform and static chemical potential which is the 
same for all electrons. Namely, the Coulomb potential 
felt by an electron from other electrons which are sepa- 
rated from it by a distance larger than Lq is described 
by a part of Eq. (7) given by 



^ roc f-OO 

= -/ dxp(x)p [ dyU(x-y) 

* J — oo J x-\-Lq 

+ f X L ° dyU(x-y)}, (28) 

J — oo 



where po is the average density. Eq. (28) represents a 
uniform one-body potential. On the other hand, when 
two electrons are closer to each other than Lq that their 
positions are correlated; then their Coulomb interaction 
has to be included in the quadratic part of the bosonic 
Hamiltonian (Eq. (10) to be specific) which governs the 
density fluctuations. In short, portions of the the sys- 
tem which are separated by distances larger than Lq are 
uncorrelated, whereas, within a distance Lq, the system 
can be described in terms of a TLL with K « 1. 

Before ending this section, we would like to present 
a simple physical understanding of the important result 
in Eq. (26) which does not make use of the replica 
idea. First, let us consider a weak ^-function impu- 
rity of strength Vq at one point in a TLL. According 
to an RG equation derived by Kane and Fisher [19], 
the impurity strength V flows according to the equa- 
tion dV/dl = (1 — K)V; hence, at a distance scale L, 
we have V(L) ~ Vo (L/a) 1 ~ K . From quantum mechan- 
ics, we know that the reflection amplitude for an electron 
scattering from a single impurity is proportional to V if 



V/ vp is small. Now suppose that each site in the lattice 
has a random impurity of strength V(x) which satisfies 
< V(x)V(y) >= D^5(x — y). This means that at each 
site, V(x) is of the order of \JD^, and its sign is ran- 
dom. Over a length L, there are L/a reflections since 
there is a reflection at each site. In order that these 
L/a reflections should add up to a total reflection am- 
plitude of order 1 (which would localize the electron), 
we must have \]L/a V(L) <~ 1 (assuming that L/a is 
large, and using the well-known result for the sum of 
a large number of random terms in a Brownian mo- 
tion). This gives us the estimate L /a ~ 1/(V(L )) 2 , 
i.e., L /a ~ l/(V (L /a) 1 - K ) 2 . Since V ~ y/I%, we 
obtain L Q /a ~ (1/£> C )V(3-2J0]. 

Our approach shows that the two ways of introducing 
disorder, (i) at random positions with a density rij but 
with equal strengths Vq as in Rcf. [13], and (ii) at all 
sites but with random strengths £(#) as in the present 
case, have the same effect. This is because the only thing 
which matters finally is the disorder parameter which 
appears in the expression in Eq. (26) for the localization 
length [17]. The exact nature of the randomness is not 
relevant for determining the form of Lq. The relation 
between the two sets of disorder parameters is given by 
~ V 2 rii. Using this relation, one can write the expres- 
sion for L given in Ref. [13] directly in terms of and 
check that Lq is proportional to (l/Z^) 1 / 3 (except for a 
logarithmic correction arising due to the LRCI [14]). This 
is also what the RG equation yields for 77 >> 1. Hence, 
our approach provides an independent (RG) derivation 
of Lq which is an alternative to the derivation using con- 
siderations of energy minimization [14,15]. We have also 
presented above an alternative understanding of the ex- 
pression in Eq. (26) which does not use the replica idea 
and is based only on the RG equation for a single impu- 
rity. 



III. DENSITY OF STATES 

We are now ready to discuss the density of states. For 
energies satisfying uj > vp/L , the disorder parameter 
D is small, and our system can be described by a clean 
TLL. The DOS of a TLL is given by the relation [9] 



D(lu) ~ 



(1 - K) 2 = VTTy 
2K ~ 2 



+ 



1 



(29) 



where we have used the relation K = l/\/l + rj. This is 
the result quoted in Eq. (4). However, we now have to 
determine what value of L one should take in Eq. (19) in 
order to determine the value of rj to be used in Eq. (29). 

It turns out that the value of u> itself determines the 
appropriate value of rj to use in Eq. (29). If u) = Vf/Lq, 
then the length scale of interest is vp/to which is equal 



4 



to Lq. The value of 77 to use in Eq. (29) is then given 
self-consistently by Eqs. (19) (with L replaced by L ), 
(21) and (26). However, if lo > vf/Lq, then the length 
scale of interest is L — vp/uj which is smaller than Lq. 
Then the Coulomb interaction should be cut-off at the 
length scale L as indicated in Eq. (18). The reason for 
this lies in the basic idea behind the RG method, namely, 
that the properties of a system at a length scale L (or, 
equivalently, at an energy scale vf/L) are governed essen- 
tially by the modes whose wavelengths are smaller than 
L. (The modes whose wavelengths are larger than L only 
contribute to a uniform and static chemical potential in 
the sense described in Eq. (28)). If the Coulomb inter- 
action is cut-off at the distance L, the value of rj to be 
used in Eq. (29) is simply given by Eq. (19). We thus 
conclude that as long as u> > Vf/Lq, the DOS of states 
is given by Eq. (29) where 77 depends on u> through Eq. 
(19) with L = vp/ui. 

We note that the expression in Eq. (4) agrees with 
the semiclassical result in Eq. (3) in the limit r\ » 1. 
Moreover, Eq. (4) correctly reproduces the expression of 
the DOS of a nonintcracting disordered electronic system 
(for which 77 = and K = 1), where there is a finite DOS 
at the Fermi energy [12]. 

We now consider the case of low energies satisfying 
u> < vp/Lo- In this case, the RG equations in (25) imply 
that the disorder is strong; hence the system cannot be 
described by a clean TLL. However, this is precisely the 
regime described by the localized limit discussed in Rcfs. 
[3-5] , where the localization length Lq is smaller than the 
length scale of interest, namely, L = vp/uj. We therefore 
expect the DOS in this regime to be described by Eq. 
(2). 

We have thus found expressions for the DOS in the 
two regimes lo > vf/Lq and u> < vp/Lo (strictly speak- 
ing, lo » vf/Lq and lo « vp/L respectively). The 
DOS must of course cross over smoothly from one regime 
to the other, but this cross over is difficult to determine 
analytically. 

Let us summarize our key results. In the limit of weak 
disorder, the system is power-law correlated only over a 
large distance of order Lq, called the localization or pin- 
ning length. At that length scale, the system is described 
by a TLL in which the LRCI drastically reduces the inter- 
action parameter K to a value much smaller than unity; 
such a value of K is quite uncommon in the literature. 
(Thus the Coulomb interaction has the same effect as a 
very strong short-range interaction) ! We note that in the 
absence of disorder, a system with a LRCI is not equiv- 
alent to one with a short-range interaction. For a clean 
system, the LRCI leads to a Wigner crystal [10] for which 
a TLL description is not valid. In our formalism, we can 
see this by noting that if the disorder was absent (i.e., 
= 0), Lq would be infinite; then rj and K would re- 
spectively be 00 and zero, and the TLL description would 
no longer be valid. On the other hand, for the disordered 



but noninteracting system (i.e., q 2 = 0), we have 77 = 
and K — 1 even if Lq is very large. 

The TLL description enables us to derive the DOS eas- 
ily using well-known relations without having to make use 
of a semiclassical approach [13]. Even though the final 
expressions for the DOS in Eqs. (3) and (4) match in 
the limit of 77 » 1, the expression in Eq. (3) derived 
in Ref. [13] shows a gap at the Fermi energy of the form 
uj 1 / 2 , and thus fails to reproduce the result for the nonin- 
teracting system. Our TLL approach correctly captures 
the physics of both the interacting and noninteracting 
limits by using the appropriate values of the interaction 
parameter K. 

Finally, beyond the length scale Lq, the TLL descrip- 
tion is no longer valid since the disorder strength is of 
order 1. In this regime, the DOS is given by Eq. (2); 
thus as lo — > 0, the DOS goes to zero logarithmically 
rather than as a power-law. 



IV. DISCUSSION 

It is interesting to consider what would happen if the 
Coulomb interaction was screened. We then have 



ft-/ 



q 2 e -x/L s 



L/2 

dx , 

L/2 VX 2 +d 2 



r, = 



ttvf 



and K 



(30) 



where L s is the screening length. It is clear that Uq now 
has a finite limit, and that 77 and K no longer approach 
00 and zero respectively, as L — > 00. For the case of weak 
disorder, the localization length Lq is again given by the 
expression in Eq. (26), where K is self-consistently de- 
termined by Eq. (30) with L replaced by Lq. 

Let us consider the case in which the screening length 
L s is much smaller than Lq. Now there are three ranges 
of energy to consider for the DOS, (i) u> > vp/L s , (ii) 
Vf/L s > u> > vf/Lq, and (iii) vf/Lq > lo. In regimes 
(i) and (ii), arguments similar to those presented in Sec. 
Ill will show that the system can again be described by 
a clean TLL in which the DOS is given by Eq. (29), 
with rj being determined by lo through Eq. (30) with 
L = vf/lo. Clearly, 77 will increase as lo decreases in 
regime (i). However, 77 will stop changing appreciably 
once L becomes much bigger than L s ; hence, 77 will be 
almost constant in regime (ii). Finally, in regime (iii), 
the disorder is of order 1, the system is in the localized 
limit (the length scale of interest is larger than Lq), and 
the DOS is given by Eq. (2). 

Before ending, we would like to mention that the re- 
sults obtained here for the localization length and the 
DOS can be used to study disordered antiferromagnetic 
spin chains. In particular, the spin- 1/2 chain with ran- 
dom magnetic field and random exchanges of various 
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kinds has been studied extensively in recent years [20-22] . 
Our analysis can be applied to such systems after map- 
ping the spins to Jordan- Wigner fermions [20] . 
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